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We present preliminary results from recent high-resolution double-degenerate merger simulations 
with the Smooth Particle Hydrodynamics (SPH) technique. We put particular emphasis on verifi- 
cation and validation in our effort and show the importance of details in the initial condition setup 
for the final outcome of the simulation. We also stress the dynamical importance of including 
shocks in the simulations. These results represent a first step toward a suite of simulations that 
will shed light on the question whether double-degenerate mergers are a viable path toward type 
la supernovae. In future simulations, we will make use of the capabilities of the NuGrid collabo- 
ration in post-processing SPH particle trajectories with a complete nuclear network to follow the 
detailed nuclear reactions during the dynamic merger phase. 
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1. Introduction 

Double-Degenerate (DD) mergers are an alternative to the single-degenerate progenitor sce- 
nario for type la supernovae (SN) that has recently received a surge in attention (e.g. Saio & 
Nomoto, 2004; Yoon et al., 2007). In this process, two white dwarfs (WD) merge, pushing one 
of them over the Chandrasekhar limit, which causes it to contract and finally explode. One of the 
main advantages of this model is that the predicted rate of such DD merger events is consistent with 
the observed SN rate, in contrast to many other SNla progenitor models. This DD merger scenario 
was early on ruled out by Nomoto & Kondo (1991). They pointed out that a CO WD accreting 
at rates higher than ~ IO^^MqS^^ would undergo "accretion induced collapse", turning it into an 
Oxygen/Magnesium/Neon WD. When pushed over the Chandrasekhar limit, this WD would then 
be able to cool sufficiently through neutrino emission to avoid a type la SN, forming a neutron star 
instead. 

However, Nomoto & Kondo's study assumed constant accretion rates, and did not account 
for a complex rotation profile. Promising new results have shed a new light on their paradigm, 
opening up small windows of opportunity for SNla as a result of DD mergers (Yoon et al., 2007). 
What is now needed are realistic high-precision simulations to determine if DD mergers are indeed 
a viable path to SNla. Our Goal is to conduct a rigorous verification and validation process and 
produce high-precision SPH simulations to answer these questions, using a modified version of the 
SNSPH code (Fryer et al., 2006; Diehl et al., 2008a). In this work, we put particular emphasis on 
the importance of appropriate initial conditions and the effects of shocks in the equation of state. 

2. The Importance of Initial Conditions 

We use the self-consistent field method to setup up our initial conditions. This method was first 
developed by Hachisu (1986) and employed by Motl et al. (2002), D'Souza et al. (2006) and Motl 
et al. (2007). It iteratively solves the equilibrium configuration for unequal mass close binaries in 
co-rotation. Its biggest advantage is that initial conditions are in perfect equilibrium, imposing very 
few artifacts in the simulations due to an imprecise initial setup. Its main drawback at the moment 
is that one is bound to choose a rather simple equation of state to produce polytropes, though work 
is underway to generalize this method for a more realistic equation of state. 

In order to be able to model a rather gentle Roche lobe overflow correctly with a particle based 
method like SPH, we have developed a new SPH particle setup method that allows us to increase 
the resolution in the outer layers, while keeping an optimal distribution of particles. This new setup 
methodis based on weighted Voronoi tesselations (WVT), allows an arbitrary spatial configuration 
without imposing a lattic geometry, and minimizes particle noise and density fluctuations (Diehl 
et al., 2008b). Figure [l| shows an example setup for a <7 = 0.4 mass ratio binary. The accretor is 
set up with a constant particle density on the left, while the donor on the right (filling 99.7% of its 
Roche lobe) has an 8x higher resolution in its outer layers. 

We also note a strong difference in behavior for corotating vs. non-corotating systems. Non- 
corotating systems tend to transfer a significant fraction of the orbital angular momentum into 
spinning up stars. This results in the orbit shrinking, bringing the star into deeper Roche lobe 
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Figure 1: Initial Condition Setup for aq ^ 1.3 mass ratio white dwarf binary system. 

contact and artificially shortening the merger time scale. For more details, refer to Fryer & Diehl 
(2008). 

3. The Importance of the Equation of State: Shocks vs. No Shocks 



We have conducted various runs with different equations of state but otherwise identical initial 
conditions to determine the effect of the equation of state on the dynamics of the interaction. Figure 
^ shows snapshots of two simulations with different equations of state at the same time in the 
simulation. The left side shows a polytropic equation of state, such that the pressure is a simple 
function of the density P = kp^ . During the simulations, this essentially results in the absence of 
shocks, since the entropy of the gas is fixed on the adiabat that is described by the constants kj) and 
kA for donor and accretor. This setup is identical to that of the ^ = 1.3 mass ratio setup of D'Souza 
et al. (2006). 

The right side of Figure ^ shows an ideal gas equation of state run instead. This simulation 
now includes shocks, which dramatically changes the dynamics of the simulation. Gas accreting 
from the donor (top) is strongly shocked when it hits the surface of the accretor (bottom). The 
major difference between the two runs is that the ideal gas equation of state builds up a hot halo of 
shocked material around the accretor which quickly starts to engulf the binary, essentially forming 
a common envelope system with different dynamics. Also note how the shocked gas in the ideal 
gas run is blown off the backside of both stars. This fundamentally changes the dynamics of the 
merger, since the expelled mass carries angular momentum outward with it. The next section will 
explore more details on this shocked run. Our goal is to produce quantitative comparison with the 
grid-based simulations of D'Souza et al. (2006) and Motl et al. (2007). 



4. Case Study: A Merger with a = 1.3 Mass Ratio with Shocks 

Figure ^ shows the time evolution of a double-degenerate merger (rotating anti-clockwise) 
between two polytropes with a mass ratio of ^ = 1.3. This sequence shows the crucial 1.5 orbits 
where the most dynamic phase of the actual merger takes place (starting at around 9.5 orbits). This 
run assumes an ideal gas equation of state and includes shocks, which can be easily seen as strong 
density discontinuities in the snapshots. 
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Figure 2: Comparison of the same time step of two otherwise identical runs for an initial mass ration of 
q = 1.3, only differing in their equations of state. SPH particles are colored according to their density, 
ranging logarithmically from 10^"* to lO". The left side uses the the polytropic equation of state, which 
keeps entropy constant during the accretion process. The right side uses an ideal gas equation of state to 
include the effects of shocks. The donor (top) material now gets heavily shocked as it hits the accretor 
(bottom), which leads to the buildup of a halo that engulfes both stars. Due to our new setup method, we 
resolve the accretion stream with well over a thousand SPH particles on average. Also note that the binary 
now loses significant amounts of mass on the backside of both white dwarfs. 




Figure 3: SPH simulation of the ^ = 1.3 merger simulation with an ideal gas equation of state. The figure 
show a time sequence starting at the onset of the merging process after about 9.5 orbits. Note the importance 
of the strong shocks. The last panel shows the core of the relaxed merger remnant ('^16 orbital periods). 
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The last panel on the lower right shows the merger remnant at a much later stage. Note that the 
remnant settled into an almost spherically symmetric configuration, and is differentially rotating, 
with a fast rotating core and a hot envelope. At the end of our simulation, the extent of the halo is 
around 100 solar radii, and is still expanding. However, we do emphasize that the actual behavior 
may very well quantitatively differ when a realistic equation of state is employed. 

5. Conclusions and Outlook 

We are on the path to produce high-precision 3D SPH simulations of double-degenerate merg- 
ers with various mass ratios. We are carrying out a detailed verification (code comparison with grid 
codes, and numerical convergence study) and vaUdation (comparison to R Coronae BoreaUs stars) 
effort. Particular emphasis is put on the initial condition setup, and the choice of the equation of 
state, which we both find to strongly affect the dynamics of the merger 

Our goal is to implement more realistic equations of state to use a stellar evolution model, and 
to follow dynamically important nuclear reactions and their energy input in the simulations. As 
part of the NuGrid collaboration \ we will also use the new post-processing tool tppnp to follow 
the evolution of SPH particles with a complete nuclear network, and use this output to vaUdate our 
simulations with abundance measurements observed in R Coronae Borealis Stars (Clayton et al., 
2007). The ultimate goal of this project is to find out whether DD mergers are valid paths to type 
la supernovae and R Coronae Borealis Stars. 

References 

Clayton, G. C, Geballe, T. R., Herwig, R, Fryer, C, & Asplund, M. 2007, ApJ, 662, 1220 

Diehl, S., Fryer, C, & Herwig, F. 2008a, in Astronomical Society of the Pacific Conference Series, 
Vol. 391, Hydrogen-Deficient Stars, ed. A. Werner & T. Rauch, 221-1- 

Diehl, S., Rockefeller, G., & Fryer, C. L. 2008b, MNRAS, submitted 

D'Souza, M. C. R., Motl, P M., TohUne, J. E., & Frank, J. 2006, ApJ, 643, 381 

Fryer, C. L., & Diehl, S. 2008, in Astronomical Society of the Pacific Conference Series, Vol. 391, 
Hydrogen-Deficient Stars, ed. A. Werner & T. Rauch, 335 

Fryer, C. L., Rockefeller, G., & Warren, M. S. 2006, ApJ, 643, 292 

Hachisu, I. 1986, ApJS, 62, 461 

Motl, P M., Frank, J., TohUne, J. E., & D'Souza, M. C. R. 2007, ApJ, 670, 1314 
Motl, P. M., Tohline, J. E., & Frank, J. 2002, ApJS, 138, 121 
Nomoto, K., & Kondo, Y. 1991, ApJ, 367, L19 
Saio, H., & Nomoto, K. 2004, ApJ, 615, 444 

Yoon, S.-C, Podsiadlowski, P, & Rosswog, S. 2007, MNRAS, 380, 933 

^http : / /forum. astro . keele . ac . uk : 8080/nugrid 



5 



